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LONG-TERM  GOALS 

The  long  range  goal  of  this  collaboration  is  to  develop  a  better  understanding  of  zooplankton 
population  dynamics,  especially  at  the  boundaries  of  distributions,  through  the  use  of  specifically- 
designed  field  efforts  utilizing  bio-acoustics,  biochemical  and  molecular  analyses,  and  process 
models. 

OBJECTIVES 

We  are  examining  how  populations  of  euphausiids  maintain  high  concentrations  in  coastal  regions 
of  the  California  Current  despite  the  highly  advective  field.  We  are  testing  the  hypothesis  that 
mesoscale  eddies  and  coastal  counter  currents  create  retention  cells,  where  growth  and  reproduction 
are  rapid  and  mortality  is  reduced.  The  immediate  objective  is  to  produce  a  data  set  on  the  biomass, 
population  genetics,  and  physiological  condition  of  the  targeted  euphausiids  across  the  boundary 
between  the  eutrophic  inshore  of  the  California  current  and  the  oligotrophic  central  gyre  waters. 
This  work  is  supported  by  ONR  Biological  Oceanography. 
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APPROACH 


We  are  using  a  suite  of  techniques  -  including  molecular  genetics,  enzyme  biochemistry,  and  bio¬ 
acoustics  -  to  examine  the  mechanisms  that  define  the  boundaries  of  zooplankton  species' 
distributions  in  the  California  Current.  We  are  analyzing  the  physiological  and  population  genetic 
characteristics  of  the  target  euphausiid  species  in  samples  collected  along  CalCOFI  line  83.3.  The 
physiological  characteristics  are  enzyme  activities  and  the  genetic  characteristics  are  molecular 
variants  of  temperature-sensitive  enzymes.  The  field  work  is  accompanied  by  the  formulation  of  a 
set  of  models  of  these  animals  that  explicitly  considers  selection  on  nuclear  genetic  traits  and  the 
effect  of  this  selection  on  population  structure. 

WORK  COMPLETED 

Molecular  genetics:  Individual  Nematocelis  difficilis  were  identified  in  the  alcohol-preserved 
samples  from  the  15  MOCNESS  tows  taken  during  RR-9610,  removed  to  individual  vials,  and 
prepared  for  molecular  analysis.  A  450  base-pair  region  of  the  mitochondrial  cytochrome  oxidase 
I  (COI)  gene  was  sequenced  for  a  total  of  1 10  individuals.  The  sequence  data  have  been  analyzed 
for  evidence  of  population  genetic  structure  by  several  tests,  including  a  hierarchical  analysis  of 
molecular  variation  (AMOVA;  Excoffier  et  al.,  1992)  comparing  variation  within  and  between 
samples  and  regions  (i.e.,  coastal,  transition,  and  offshore). 

After  some  deliberation  about  technical  approaches  to  identify  a  trait  encoded  in  the  nuclear  genome 
that  is  likely  to  be  subject  to  selection,  we  decided  to  use  differential  display,  reverse  transcriptaise 
PCR  (dd-rtPCR)  to  examine  genes  that  are  differentially  expressed  (i.e.,  turned  on-or-off)  in  different 
populations  and  environmental  conditions.  We  have  extracted  usable  mRNA  fi-om  individually- 
frozen  N.  difficilis  collected  during  RR-9610,  and  are  currently  adapting  established  dd-rtPCR 
protocols  to  euphausiids.  The  dd-rt  PCR  assays  require  only  a  small  amoimt  of  tissue;  the 
biochemical  and  genetic  analyses  were  thus  done  using  the  same  individuals  (see  below). 

Biochemical  analysis:  Enzyme  kinetic  assays  were  completed  for  100  individuals  of  N.  difficilis 
collected  during  RR-9610  (see  Clark  et  al.,  1992;  Clark  and  Walsh,  1993).  The  enzymes  assayed 
were:  citrate  synthetase  (CS),  lactate  dehydrogenase  (LDH),  phosphoglucose  isomerase  (PGI),  and 
hexokinase  (HK).  In  addition,  the  same  analyses  were  done  for  300  individuals  of  several  other 
euphausiid  species.  Protein  analyses  were  conducted  on  all  samples;  all  enzyme  values  were  scaled 
to  protein.  These  data  will  be  analyzed  to  determine  the  physiological  condition  and  growth  of 
individual  euphausiids  across  the  sampled  domain. 

High-frequency  acoustics:  In  order  to  groundtruth  the  bioacoustics  data,  the  taxonomic  composition 
and  size  fi-equency  distribution  of  zooplankton  were  determined  for  MOCNESS  samples  from  RR- 
9610  by  silhouette  analysis  (see  Wiebe  et  al.,  1996).  Silhouette  analysis  was  completed  for  four  net 
samples  (spanning  the  top  100  m)  from  each  of  four  MOCNESS  tows.  The  processing  of  the 
bioacoustic  backscatter  data  continue.  The  acoustic  results  will  provide  the  biomass  context  for  the 
single-species  molecular  and  biochemical  analyses. 
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Modeling:  A  schematic  of  a  model  to  treat  the  impact  of  enzyme  expression  on  a  planktonic 
organism  is  finished.  The  process  of  simplifying  the  model  and  testing  it  for  compatibility  with  the 
observations  has  begun.  A  set  of  physical  models  covering  upwelling  systems  was  developed  based 
on  the  Miami  Community  Ocean  Model  (MICOM).  Three  types  of  plankton  models  are  currently 
up  and  running  in  these  simulations.  The  first  two  are  a  simple  NPZD  model  and  a  multiple 
foodchain  model  that  is  also  constructed  with  nutrients  as  a  conservative  currency.  These  models 
are  used  as  a  background  age  (stage)  and  metabolically  structured  model  that  runs  on  particles,  and 
can  also  simulate  genetic  characteristics  of  planktonic  populations. 

RESULTS 

Oceanographic  setting:  The  hydrographic  analysis  from  the  domain  sampled  during  RR-9610 
indicated  the  presence  of  a  coastal  eddy  transected  by  CalCOFI  line  83.3,  which  was  the  premise  of 
our  proposal.  Based  on  surface  hydrography  and  circulation,  our  transect  of  15  MOCNESS  tows 
sampled  from  offshore,  nutrient-poor  waters;  a  transition  region  of  intermediate  properties;  and 
coastal,  nutrient-rich  waters  (see  RR-9610  cruise  report,  T.  Hayward,  ed.). 

Molecular  genetics:  We  have  characterized  the  molecular  genetic  diversity  and  population 
structure  of  Nematoscelis  difficilis  across  the  sampled  domain  using  DNA  sequence  variation  of  the 
mitochondrial  gene,  cytochrome  oxidase  I  (COI;  see  Bucklin  et  al.,  1997).  Interestingly,  the  species 
population  is  heterogeneous  across  the  California  Current  (Fig.  1).  The  sample-to-sample 
heterogeneity  (P  <  0.015  based  on  a 
chi-square  test  of  haplotype 
frequencies)  suggests  a  complex 
response  to  the  environmental  gradient 
along  the  transect.  Based  on  AMOVA, 

8%  (P  <  0.032)  of  the  variation  was 
explained  by  variation  between 
samples;  there  was  no  evidence  of 
regional  groupings  (i.e.,  onshore  vs 
offshore)  in  the  spatial  patterns  of 
mtDNA  diversity. 

We  have  modified  protocols  for  assay 
of  gene  expression  differences  by 
differential  display-reverse 
transcriptase  PCR  (dd-rtPCR)  for 
euphausiids,  and  have  obtained  high-  Figure  1.  MtCOI  haplotype  frequencies  of  the  euphausiid 
quality  mRNA  from  individual  N.  in  samples  collected  along  CalCOFI  line  83.3; 

difficilis.  We  are  now  screening  for  numbers  are  CalCOFI  station  designations.  Haplotype  frequencies 
gene  expression  differences  between  differed  significantly  among  the  samples  (P<0.015),  but  the  spatial 
individuals  collected  from  the  coastal  heterogeneity  of  the  populations  was  not  resolved  into  regional 
eddy  vs.  those  collected  offshore  using  patterns. 


differential  display.  In  addition,  we  are  working  to  obtain  cDNA  sequences  for  glycolytic  enzymes, 
including  phosphoglucose  isomerase  (PGI)  and  others  assayed  biochemically.  Comparison  of  the 
DNA  sequence  of  the  encoding  gene  with  the  assays  of  enzyme  activities  for  the  same  individuals 
will  allow  us  to  determine  whether  any  observed  differences  in  individuals  collected  along  the 
environmental  gradient  result  from  molecular  genetic  diversity  (i.e.,  DNA  sequence  variation)  or 
physiological  adaptation  (i.e.,  differences  in  gene  expression  and  enzyme  concentrations). 

Biochemical  analysis:  Frozen  samples  of  N.  difficilis  yielded  high-quality  data  on  the  physiological 
condition  of  individual  euphausiids,  based  on  the  kinetic  activities  of  citrate  synthetase  (CS), 
hexokinase  (HK),  lactate  dehydrogenase  (LDH),  and  phosphoglucose  isomerase  (PGI).  Activities 
of  all  enzymes  were  measurable,  high,  and  highly  variable.  Although  enzyme  concentrations  were 
not  significantly  different  among  individuals  collected  in  different  MOCNESS  tows  (by  analysis  of 
variance:  CS,  P<0.54;  HK,  P<0.14;  LDH,  P<0.25;  PGI,  P<0.76),  it  is  possible  that  patterns  of 
enzyme  concentration  may  be  related  to  genetic  characters.  We  will  expand  our  statistical  analyses 
of  the  molecular  and  biochemical  data  to  include  multivariate  analysis  of  DNA  sequence  variation 
of  mtCOI  and  enzyme-encoding  genes 
and  differential  gene  expression,  once 

all  the  analyses  are  complete.  ^2.4  Stn  77.4 

High-frequency  acoustics:  Silhouette 
analysis  of  taxon-specific  biomass 
revealed  considerable  variation  in  the 
taxonomic  composition  of  samples 
along  the  transect  (Fig.  2),  suggestive  of 
mesoscale  patchiness  structure.  The 
silhouette  analysis  data  will  provide  the 
biomass  context  for  our  single-species 
modeling  and  will  be  used  to  ground- 
truth  the  bioacoustical  observations. 

Modeling:  The  numerical  model 

incorporating  population  genetic  data 
has  been  developed  (Olson  et  al.,  1998) 
and  is  being  modified  for  the  California 
Current  coastal  region.  The  mtDNA  Figure  2.  Biomass  fractions  of  copepods,  euphausiids,  and  salps 
haplotypes  are  used  a  tags  of  dispersal  determined  from  silhouette  analysis  for  depth-stratified  samples  from 
in  the  numerical  model;  previous  efforts  four  MOCNESS  tows  at  standard  stations  (Stn)  along  CalCOFI  line 
have  demonstrated  the  model  s  ability  to  83.3.  All  other  taxa  are  grouped  as  "other".  These  data  are  important 
explain  and  reproduce  complex  to  understand  the  biomass  contribution  of  euphausiids  and  to  ground- 
geographic  patterns  of  mtDNA  truth  the  bioacoustical  analysis.  Data  for  the  0-25m  sample  at  Stn 
haplotype  frequencies  (Olson  et  al.,  77.4  are  not  available  (N/A). 

1996,  1998). 
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A  preliminary  study,  based  on  transport  of  a  planktonic  copepod  around  the  N.  Atlantic  gyre, 
revealed  significant  spatial  pattern  in  the  population  genetic  character  of  zooplankton  as  they  are 
advected  through  varying  environmental  carrying  capacities.  The  model  results  were  consistent  with 
the  observations  in  the  Gulf  Stream  system,  based  on  a  published  study  of  the  copepoda, 
Nannocalanus  minor  (Bucklin  et  al.,  1996).  A  manuscript  submitted  to  the  Journal  of  Marine 
Research  is  now  being  revised  for  resubmission. 

Synthesis:  We  are  working  toward  synthesis  of  our  diverse  data  sets,  and  full  assimilation  into  the 
numerical  population  dynamic  models.  We  continue  to  refine  our  integrated  tools  (genetics, 
bioacoustics,  and  numerical  modeling)  in  order  to  evaluate  both  model  results  and  field  observations. 
We  are  also  experimenting  with  predictive  capabilities  in  various  field  regions,  including  the  coastal 
boundary  region  of  the  California  Current. 

IMPACT 

Our  fundamental  approach,  which  combines  experimental  data  and  modeling,  should  provide  new 
insights  into  interpretation  of  population  dynamics  in  the  ocean.  In  addition,  the  integrated  analysis 
of  biochemical,  molecular,  and  high-frequency  acoustic  data  will  provide  new  information  for 
biological  oceanographers  who  seek  to  understand  complex  processes  and  will  also  provide  new 
indices  for  rapid  assessment  of  the  planktonic  assemblage.  The  impact  of  this  work  will  be  to 
improve  biological  models  of  ocean  assessment  and  prediction. 

TRANSITIONS 

The  model  codes  will  be  made  widely  available  after  they  have  been  further  tested  and  reviewed  for 
publication.  Olson  is  discussing  possible  transitions  of  the  Lagrangian  code  to  Naval  Research  Lab 
(NRL)  via  collaboration  with  Dr.  John  Kindle  at  NRL. 

RELATED  PROJECTS 

Using  samples  collected  during  RR-9610,  A.B.  has  examined  population  genetic  diversity  and 
structure  of  the  copepoda,  Calanus  pacificus,  in  coastal  regions  of  the  California  Current.  These  data 
are  being  analyzed  now,  to  determine  whether  the  nearshore  regions  provide  opportunities  for 
retention  and  persistence  for  one  of  the  most  abundant  and  ecologically  important  members  of  the 
California  Current  zooplankton  assemblage.  This  analysis  will  help  us  interpret  our  results  for 
Nematoscelis  difficilis,  which  has  different  behavior  patterns  and  preferred  habitats,  but  may  also  be 
retained  in  coastal  cells  and  counter  currents. 
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